NEW FEATURES OF THE PHASE TRANSITION TO 
SUPERCONDUCTING STATE IN THIN FILMS 



D. V. SHOPOVA*, T. P. TODOROV* 

CPCM Laboratory, Institute of Solid State Physics, 
Bulgarian Academy of Sciences, BG-1784 Sofia, Bulgaria. 

* Corresponding author: sho@issp.bas.bg 

' Permanent address: Joint Technical College at the Technical University of Sofia. 

Key words: fluctuations, latent heat, order parameter, equation of state. 
PACS: 05.70.Jk, 74.20.De, 74.78.Db, 

Abstract 

The Halperin-Lubensky-Ma (HLM) effect of a fluctuation-induced change of the 
order of phase transition in thin films of type I superconductors with relatively small 
Ginzburg-Landau number k is considered. Numerical data for the free energy, the 
order parameter jump, the latent heat, and the specific heat of W, Al and In are 
presented to reveal the influence of film thickness and material parameters on the 
properties of the phase transition. We demonstrate for the first time that in contrast 
to the usual notion the HLM effect occurs in the most distinct way in superconducting 
films with high critical magnetic field H c q rather than in materials with small k. The 
possibility for an experimental observation of the fluctuation change of the order of 
superconducting phase transition in superconducting films is discussed. 

1. Introduction 

Our paper is intended to clarify the best conditions for an experimental observation of 
the Halperin-Lubensky-Ma effect (HLM) |2| of a fluctuation-induced first-order phase 
transition from normal to Meissner phase in a zero external magnetic field for thin films of 
type I superconductors [HI IH |S] . For this purpose we present new theoretical results about 
the thermodynamics of the phase transition from normal to superconducting state in a zero 
external magnetic field. 

The HLM effect is predicted theoretically jT] to occur in pure IH El and disordered [7J 
[HIE] bulk, (three-dimensional - 3D), and 2D [TOJ superconductors, as well as in quasi-2D 
superconducting films [31 Ej but up to now it has not been observed in experiments. The 
calculated effect is very small in 3D superconductors and is not possible to be detected even 
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for a high purity of the sample and perfection of the crystal lattice |U 121 HI • Recently, it 
has been shown [HI E] that in thin (quasi-2D) films the HLM effect is much stronger than 
in 3D samples and could be observed by available experimental techniques if the type of 
superconductor and film thickness are properly chosen for the experiments; for a review, 
see also, Refs. [HJH2]. This result gives an opportunity to search for the effect in suitable 
superconducting films. 

The HLM effect appears as a result of the interaction between the superconducting order 
parameter if)(x) and the vector potential A(x) of the magnetic induction in the Ginzburg- 
Landau (GL) free energy of a superconductor. According to the theoretical paradigm intro- 
duced for the first time in the scalar electrodynamics by Coleman and Weinberg (CW) [T3*j . 
this effect should occur in all physical systems described by Abelian-Higgs models where 
a scalar gauge field (like ip in superconductors) interacts in a gauge invariant way with 
another vector gauge field (like the vector potential A in superconductors). In addition 
to the mentioned examples of superconductors and scalar electrodynamics, the same type 
of interaction plays an important role in the nematic-smectic A phase transition in liquid 
crystals El] and phase transitions in the early universe [Tj\. HLM effect may be 

also relevant to quantum phase transitions in superconductors fHl E3 120] and itinerant 
ferromagnets |21j . 

On the other side, there are certain theoretical investigations, based on Monte Carlo simula- 
tions [22] and the so-called "dual model" [2B] which do not confirm the fluctuation-change 
of the order of the phase transition (see also Refs. [EI). That is why, extensive experi- 
ments intended to verify the existence of the effect were made in liquid crystals; see, e.g., 
Refs. [12]. But in liquid crystals the weakly- first order phase transition predicted by CW 
and HLM can be obscured by similar effects due to the strong crystal anisotropy, while the 
recent result [3] about the considerable enhancement of HLM effect in suitable supercon- 
ducting films can be used for more reliable experiments. For this aim we need to find the 
best material parameters and the most suitable film thickness, having in mind some purely 
experimental problems that may appear. 

Recently, we partly solved the problem for Al films [3] nevertheless additional theoretical 
investigations should be done. In this paper we shall present some new theoretical predic- 
tions for Al films as well as new numerical data for thin films of W and In. The choice 
of this element superconductors is made for their relatively small GL number k = (A/£), 
which allows a more distinct appearance of the HLM effect in both bulk and thin film 
superconductors [TJ El IS] ; here A is the London penetration depth and £ is the coherence 
length [23| . 

We focus our attention on numerical data for the behavior of the free energy and directly 
measurable thermodynamic quantities like the order parameter jump, the latent heat, and 
the specific heat. A surprising result of our analysis of the data for W, Al, and In is that 
the HLM effect in thin films is stronger in case of relatively high zero-temperature critical 
magnetic field H c q rather than for relatively small GL number k, as claimed in preceding 
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papers [JJ H El IH El • 

Our investigation is based on the theoretical results from preceding papers PflEHHEJ- In 
Sec. 2 we shall outline the theoretical framework of our study. In Sec. 3 we present our 
analysis of thin films of tungsten (W), aluminium (Al), and indium (In) and a discussion 
of the results with a special emphasis on their application to experiments. In Sec. 4 we 
summarize our findings. 

2. Theoretical basis 

Our investigation is based on the Ginzburg-Landau free energy [21] of a D-dimensional 
superconductor with volume V = (Li...Ld) given by 
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where a = a (T — T c0 ) and 6 > 0, are the Landau parameters, e = |e| is the electron charge, 
ip(x) is the order parameter, and A(x) is the vector potential of the magnetic field. 

The critical temperature T c q corresponds to the second order phase transition which occurs 
in a zero external magnetic field when the fluctuations, 5ip(x) and 5A(x), of both fields ip(x) 
and A(x) are neglected. Usually, this case is considered in the low-temperature supercon- 
ductors, where the effect of the superconducting fluctuations 5<p(x) on the thermodynamics 
is very small and practically uninteresting, and the same has been supposed for the magnetic 
fluctuations SA(x) before the appearance of HLM paper [Tj. In our study the supercon- 
ducting fluctuations are ignored as negligibly small which is a suitable approximation in 
type I superconductors where A <C £. But we take into account the magnetic fluctuations 
to a full extent. Then the normal-to-superconducting phase transition in a zero external 
magnetic field turns out of first order at an equilibrium phase transition temperature T eq 
that is different from T c0 . So, our task will be to point the type of superconductors, where 
this picture may be valid and investigate the properties of the first order phase transition. 

We shall follow the theoretical approach described in details in preceding papers [TJ El El E] 
where an effective free energy F e ff(ip) of the type I superconducting film was obtained 0. 
There we also neglected the superconducting fluctuations and did the calculation of the 
mean-field value of the uniform (x- independent) superconducting order parameter ip in a 
self-consistent way after taking into account the magnetic fluctuations through an exact 
integration out of the field A(x) in the partition function of the superconductor. As the 
external magnetic field is equal to zero the regular part Aq = (A — 5 A) of A related to it 
can be set equal to zero, too. Then A = 5A, therefore, we consider the net effect of the 
magnetic fluctuations. 

Having in mind that in type I superconductors a stable vortex phase cannot occur, we 
again can assume that the order parameter ip that describes the uniform Meissner phase 
in the bulk of the superconducting film is x-independent. Our investigation is based on 
the quasi- macroscopic GL theory so we must consider films of thickness Lq ^> do, where 
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do is the lattice constant. In such films the surface energy can be ignored and one can use 
periodic boundary conditions without a substantial departure from the real situation. In 
this way, the surface effects as a source of a spatial dependence of the order parameter ip 
are also eliminated. 

Following [TJ |21 El IH E] we present the effective free energy density f(ip) = F e ^{ip)/V of a 
3D superconducting slab of volume V = (LiL 2 L Q ) and thickness L Q in the form: 

f(<p) = {2t ^ 2 + ^ + C(l + to) [(1 + W 2 ) In (1 + W 2 ) - /Vln (/V)] } , (2) 
where 

_ 2^k B T c0 

U ~ T C2U2 ■ K 6 ) 

Here <p = (\i/j\/\i/jo\) is the dimensionless order parameter defined with the help of the zero- 
temperature value \if) \ = \^{T = 0)| = (a T c0 /b) 1/2 of \if)\, t = (T-T c0 )/T c0 , \i = (£o/ti\ ) 2 
is given by the zero-temperature value £o = (^ 2 /4maoT c0 ) 1//2 of £ and Ao = (fr/po^o^co) 1 ^ 2 
is the zero-temperature penetration depth; (p = 8ire 2 /mc 2 ). We also use the notations: 
A(T) = A /|t | 1/2 and f (T) = £/|t | 1/2 - The critical magnetic field at T = is given by [21] 
H c o = a T c0 (47r/6) 1//2 . The relations of H c0 and £ with b and a , respectively, can be used 
together with the experimental data for H c0 and £o in concrete superconducting substances 
in order to calculate the parameters b and ao- 

The equilibrium order parameter <po > corresponding to the Meissner phase can be easily 
obtained from the equation df(ip)/dip = and Eq. (j2J): 
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The logarithmic divergence in Eq. (@J) has no chance to occur because (po is always positive 
and does not tend to zero. 

We shall use the notations from Ref. [S] for the entropy jump 5 s and the specific heat 
jump 5C at the equilibrium phase transition point T eq of first order corresponding to a zero 
external magnetic field. Here we shall give the previously calculated results for the leading 
terms in these quantities (terms of higher order are neglected as small), namely, 



and 



The latent heat of the phase transition is given by Q = T eq 5s and Eq. (0]). Since the 
temperatures T eq and T c0 have very close values, the difference between the values of Q, 
5s, and 5C at T c0 and T eq , respectively, can also be ignored, for example, \5C(T eq ) — 
5C(T c0 )\/5C(T c0 ) <C 1 and we can use either 5C(T c0 ) or 5C(T eq ) Here the jumps 5s, 
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8C, and Q are all taken at the equilibrium phase transition value T eq but we shall not 
supply them with the subscript "eq" as we do for other quantities. 



for negligibly small ^-fluctuations which means that they are applicable for low-temperature 
nonmagnetic superconductors (T c0 < 20K). Because in experiments the external magnetic 



T c = T C (H) < T c0 in type II superconducting films and this will obscure the HLM effect. 
Note, that the magnetic field H generates additional entropy jump at the phase transition 
point T C (H) and this effect can hardly be separated from the entropy jump (J5J) caused by 
the magnetic fluctuations in the close vicinity of T c q. Therefore, in experiments intended to 
a search of the HLM effect we must choose type I superconductors. The second important 
point is connected with the value of the square ip 2 q which is proportional to the supercon- 
ducting current (j s ~ ~ V^eqElI) an d to the equilibrium jumps 5s and 5Q. In 3D 
superconductors the ratio (Q/5C) depends on H^ ~ e c K~ 6 , where e c ~ 1CT 16 denotes the 
extremely small Ginzburg-Levanyuk critical region [llj of low-temperature superconduc- 
tors pQ. Therefore, the latent heat in these 3D superconductors can hardly be observed 
in experiments. But in thin films the substantial dependence of the entropy 5s and the 
specific heat 5C is on the critical magnetic field H 2 Q , as shown by Eqs. © - © and the 
analysis in Sec. 3. 

The equations (J2J and (jlj corresponding to quasi-2D films are quite different from the 
respective equations [TJ [2] for bulk (3D-) superconductors but it is easily seen that the 
relatively large values of the order parameter jump (p 2 in thin films again correspond to 
relatively small values of the GL parameter k. That is why we consider element supercon- 
ductors with small values of k and study the effect of this parameter, the critical magnetic 
field H c0 and the film thickness L on the properties of the fluctuation-induced first order 
phase transition. 

Theoretical results we have used in this Section for quasi-2D superconducting films are 
consistent with the theory [25] of 2D-3D crossover phenomena near phase transition points 
and the 2D-3D crossover theory [231211 of the HLM effect; see also Ref. [21]. 

3. Results and discussion 

We use experimental data for T c0 , H c0 , £ and k for W, Al, and In published in Ref. [22] (see 
Table 1). In some cases the GL parameter k can be calculated with the help of the relation 
K = (Aq/£o) an d the available data for £ an d A . In other cases it is more convenient to 
use the following representation of the zero-temperature penetration depth: 



Eqs. (j2j) - (JHJ) are valid for thin films (ao «C L$ ~ £o) in a zero external magnetic field H and 



field cannot be completely eliminated, vortex states may occur for H 
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The value of \i/jq\ in Table 1 is found from 
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Eqs. (J7J) and (jSJ) are obtained from the formulae given after Eq. (j2J). Besides we calculate 
the parameter C = CLq with the help of Eq. (J3J) and the data in Table 1. 

Note, that the experimental data vary within 5-10% depending on the experimental tech- 
nique used in measurements . Moreover, these data correspond to bulk samples and may 
differ within 10-20% from those for very thin films (L < 10 -2 /Lim). However, these varia- 
tions in the experimental data do not essentially affect our results. 



Table 1. Values of T c o, 


Hco, Co, K , 


and -00 


for W, 


Al, In. 


substance T c0 (K) 


H c0 (Oe) 


£o (M m ) 


K 


IV'ol x lO" 11 


W 0.015 


1.15 


37 


0.001 


0.69 


Al 1.19 


99.00 


1.16 


0.010 


2.55 


In 3.40 


281.5 


0.44 


0.145 


2.0 



Table 2. Values of t eq , <^ eq , and Q (erg/cm 3 ) for films of W, Al, and In with different 



thicknesses Lq (fim). 







Al 






In 






w 








teq 


^cq 


Q 


teq 




Q 


teq 


Veq 


Q 




0.05 


-0.00230 


0.041 


-1.95 


-0.00167 


0.025 


-3.94 


-0.00174 


0.039 


-1.6 x 10~ 


-4 


0.1 


-0.00147 


0.032 


-0.80 


-0.00094 


0.017 


-1.82 


-0.00118 


0.032 


-1.1 x 10- 


-4 


0.3 


-0.00070 


0.023 


-0.41 


-0.00037 


0.010 


-0.63 


-0.00064 


0.023 


- - 5.6 x HT 5 


0.5 


-0.00048 


0.016 


-0.20 


-0.00029 


0.008 


-0.40 


-0.00048 


0.020 


-4.1 x 10" 


-5 


1 


-0.00029 


0.012 


-0.11 


-0.00013 


0.006 


-0.23 


-0.00032 


0.016 


-2.7 x 10" 


-5 


2 


-0.00017 


0.009 


-0.06 


-0.00008 


0.004 


-0.10 


-0.00021 


0.013 


-1.8 x 10- 


-5 



The order parameter profile for Al films of different thicknesses is shown in Fig. 1. It is 
readily seen that the behavior of the function <po(to) corresponds to a well established phase 
transition of first order. The vertical dashed lines in Fig. 1 indicate the respective values 
of t e q — to(T eq ), at which the equilibrium phase transition occurs as well as the equilibrium 
jump <yj (T eq ) = (p eq for different thicknesses of the film. The parts of the y?o(^o)- curves 
which extend up to to > teq describe the metastable (overheated) Meissner states which can 
appear under certain experimental circumstances (see in Fig. 1 the parts of the curves on 
the r.h.s. of the dashed lines). The value of <p eq and the metastable region decrease with 
the increase of the film thickness, which shows that the first order of the phase transition 
is better pronounced in thinner films and that confirms a conclusion in Ref. |2Ti] . 

These results are confirmed by the behavior of the free energy as a function of to- We used 
Eqs. (J2J) and (|1J) for the calculation of the equilibrium free energy f[ip (t )]. The free energy 
for Al films with different thicknesses is shown in Fig. 2. The equilibrium points T eq of the 
phase transition correspond to the intersection of the f (ip )-curves with the t -a.xis. It is 
obvious from Fig. 2 that the temperature domain of overheated Meissner states decreases 
with the increase of the thickness Lq. 

The shape of the equilibrium order parameter <po in a broad vicinity of the equilibrium 
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Figure 1: Order parameter profile <p(to) of Al films of different thicknesses: L = 0.05 /xm 
("+"-line), L = 0.1 fim (o), and L = 0.3 fim (•). 
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Figure 2: The free energy /(to) f° r Al films of thickness: Lq = 0.05 fim ("+"-line), Lq 
O.f fim (o), L = 0.3 /xm (•). 
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Figure 3: Order parameter profile <p(to) of films of thickness L = 0.05 /im: W ("+"-line), 
Al (o), and In (•). 

phase transition of thin films (L = 0.05/im) of W, Al, and In was found from Eq. (0J. The 
result is shown in Fig. 3. The vertical dashed lines in Fig. 3 again indicate the respective 
values of t eq = t (T cq ), at which the equilibrium phase transition occurs as well as the 
equilibrium jump ip (T eq ) = ip eq in the different superconductors. 

The order parameter jump at the phase transition point of In (the In curve is marked by 
points in Fig. 3) is relatively smaller than for W, and Al, where the GL parameter has much 
lower values. The same is valid for the metastability domains; see the parts of the curves 
in Fig. 3 on the left of the vertical dashed lines. It is obvious from Fig. 3 and Table 2 that 
the equilibrium jump of the reduced order parameter (p eq of W has a slightly smaller value 
than that of Al although the GL number k for W has a ten times lower value compared 
with k of Al. Note, that in Fig 3 we show the jump of (p eq , but the important quantity 
is \ip\eq = l^olfeq- Using the data for L$ = 0.05/im from Tables 1 and 2 we find for iV'leq 
the following values: 0.1 x 10 11 for Al, 0.05 x 10 11 for In, and 0.02 x 10 11 for W. This 
result shows that the value of the critical filed H c q is also important and should be taken 
into account together with the smallness of GL number when the maximal values of the 
order parameter jump are looked for. Thus the value of the order parameter jump at the 
fluctuation-induced phase transition is maximal provided small values of the GL parameter 
k are combined with relatively large values of the critical field H cQ . In our case Al has the 
optimal values of these two parameters. 

The shift of the phase transition temperature t eq = |(T eq — T c0 )|/T c0 , the reduced value 
ip cq of the equilibrium order parameter jump \ip\ e q, and the latent heat Q of the equilib- 
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rium transition are given for films of different thicknesses and substances in Table 2. The 
thicknesses are chosen so as to ensure the validity of the theory (HI ED] used in our analysis 
and to satisfy other important requirements presented in Sec. 4. The data in Table 2 show 
that the shift of the phase transition temperature is very small and can be neglected in all 
calculations and experiments based on them. The values for ip eq for different L and those 
for | -00 1 given in Table 1 confirm the conclusion which we have made for films of Al, In, and 
W with Lq = 0.05/im. The latent heat Q has maximal values for In, where the critical field 
is the highest for the considered materials. 

4. Conclusion 

In contrast to our initial expectations that films made of superconductors with extremely 
small GL parameter k such as Al and, in particular, W will be the best candidates for 
an experimental search of the HLM effect, our careful analysis definitely gives somewhat 
different answer. The Al films still remain a good candidate for transport experiments 
through which the jump of the order parameter at the phase transition point could be 
measured but surprisingly the W films turn out inconvenient for the same reason because 
of their very low critical field H c q. Although In has ten times higher GL number k than 
Al, the In films can be used on an equal footing with the Al films in experiments intended 
to prove the order parameter jump. Here the choice of one of these materials may depend 
on other features of experimental convenience. As far as caloric experiments are concerned, 
the In films seem the best candidate for their high latent heat. 

We have presented the theoretical justification and predictions intended to support exper- 
iments on the observation of magnetic fluctuations and HLM effect near the normal-to 
superconducting transition in a zero external magnetic field. Besides, we have demon- 
strated for the first time that the experiments can be most successfully performed in type 
I superconductors with relatively high critical magnetic fields. 
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